Spin-charge separation in ultra-cold quantum gases 
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We investigate the physical properties of quasi- ID quantum gases of fermion atoms confined 
in harmonic traps. Using the fact that for a homogeneous gas, the low energy properties are 
exactly described by a Luttinger model, we analyze the nature and manifestations of the spin- 
charge separation. Finally we discuss the necessary physical conditions and experimental limitations 
confronting possible experimental implementations. 



One dimensional (ID) quantum liquids are very rich 
and interesting systems. In spite of their apparent con- 
ceptual simplicity, both the ground state and the excita- 
tions exhibit strong correlation effects and posses a num- 
ber of exotic properties, ranging from spin-charge separa- 
tion to fractional statistics (see 0, ^ and ref. therein). 
Progress in creating, manipulating and studying ultra- 
cold quantum gases with controlled and adjustable in- 
teractions P|, and in particular the recent develop- 
ment of ID magnetic and optical wave guides opens the 
door for a new and clean physical realization of such ID 
systems with the tools of atomic physics and quantum 
optics. While most of the recent theoretical and experi- 
mental work has focused on ID Bose gases (as a Tonks 
gas or a quasicondensate) |^ progress in cooling Fermi 
gases into the quantum degenerate regime [Q point to 
the possibility of realizing a Luttinger Liquid (LL) 
with cold fermionic atoms. One of the key predictions 
of the Tomonaga-Luttinger model (LM) for interacting 
fermions is spin-charge separation 0]. It is a feature 
of interacting spin- 1/2 particles and manifests itself in 
complete separation in the dynamics of spin and den- 
sity waves. Both branches of the excitations are sound- 
like and characterized by different propagation velocities. 
This phenomena is one of hallmarks of a Luttinger liq- 
uid, however it has never been been demonstrated in a 
clean way in an actual condensed matter system (see e.g. 
0, |l|). It is the purpose of this Letter to analyze in 
detail the conditions of realizing an (inhomogeneous) LL 
with a gas of cold fermionic atoms in ID harmonic trap 
geometries, and in particular to study the possibilities 
of seeing spin-charge separation in the spectroscopy and 
wave packet dynamics of laser excited ID Fermi gases. 

The simplest example of a Luttinger liquid made of 
a gas of cold atoms consists of fermionic atoms with 
two ground states representing a spin-1/2 under quasi-lD 
trapping conditions. We assume the atoms to be cooled 
below the Fermi-degeneracy temperature ksTp ~ Nhw, 
where N is the number of particles and uj is the frequency 
of the longitudinal confinement. The condition for a 
quasi- ID system is tight transverse trapping in an exter- 
nal potential with the frequency lo_\^ exceeding the char- 
acteristic energy scale of the longitudinal motion. Due 
to the quantum degeneracy the longitudinal motion has 
all the energy levels up to the Fermi-energy ep ^ ksTp 
filled. Thus we require the total number of particles to 



Fig. 1. (a) Schematic setup: a two component fermi gas is 
trapped in a harmomic potential in a ID configuration. At 
time t = a short laser pulse focused near the center of the 
trap excites (i) a density (charge) or (ii) a spin-wavepacket. 
(b) Wave packet dynamics for different times as function of 

position (in units of the Thomas Fermi radius Rtf): 
spin-charge separation manifests itself in a spatial separation 
of the spin (solid line) and density (dashed line) wave 
packets (shown at half a trap oscillation period cut = 0.5), 
which can be probed by a second short laser pulse at a later 
time. The parameters correspond to A = lO'^ ^Li atoms in a 
trap with ui = IHz (with coupling parameter ^ = 1, see text). 



be restricted by A^ ^ uj±/u!, which for realistic traps 
is of the order of a few hundred or thousand. Let us 
now turn to an estimate of the conditions to reach the 
strongly interacting limit. Since the interaction between 
the atoms has a range much smaller than the interparti- 
cle spacing, at low temperatures only collisions between 
particles with different spins are allowed by the exclu- 
sion principle. Therefore, all the relevant interactions 
are characterized by a single parameter, the scattering 
length a corresponding to inter-component interaction. 
The effective ID interaction can thus be represented as 
a zero-range potential of the strength g — ^tx^t} a j ml\^ 
where is the width of the ground state in the trans- 
verse direction (a ^ ^j_)|^ and m is the mass of the 
gas particle. The interaction strength in a Luttinger liq- 
uid is then characterized by the dimensionless parameter 
^ = g/irhvp, where mvp/2 = kBTp is the Fermi ve- 
locity (at trap center). Remarkably, in a trapped gas 
^ ~ {a/l±){uj±/LuNy/^ and thus can be tuned exter- 
nally either by changing the transverse confinement, or 
by changing the scattering length by magnetic field. The 
ratio of the transverse and the longitudinal frequencies is 
quite large, we can easily reach the strong coupling limit 
^ ~ 1, even in a dilute gas. 

Below we will study spin-charge separation according 
to the schematic setup outlined in Fig. 1. We assume that 
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a short far ofF-resonant laser pulse is focused at the center 
of the harmonic trap with a two-component atomic Lut- 
tinger liquid, where depending on the laser parameters 
(e.g. light polarization) density and spin wave packets 
can be excited (Fig. la). Spin-charge separation mani- 
fests itself in different propagation speed of the spin and 
density wave packets (see Fig. lb). This can be probed 
at a later time with a second short laser pulse. Spin de- 
pendent optical potentials can be generated by a laser 
tuned e.g. between fine structure levels of excited Alkali 
states acts on the ground state "spin" in a way equiva- 
lent to the external magnetic field interacting with the 
spin density and thus introduce a spin density perturba- 
tion. The goal of the following derivations is thus to (i) 
derive the frequencies of the spin and charge modes of 
atomic Luttinger liquid confined in a harmonic trap, and 
(ii) to discuss the wavepacket dynamics as superposition 
of these modes. 

We analyze the properties of a trapped Luttinger liquid 
by combining Haldane's low energy hydrodynamic de- 
scription Q with a local density approximation. As was 
shown by Haldane, any homogeneous interacting quan- 
tum liquid in one dimension can be described by an af- 
fective hydrodynamic Hamiltonian, which completely de- 
scribes the behavior at wavelengths much larger than the 
interparticle spacing. In the case of spin-1/2 fermions, 
the effective Hamiltonian reads 
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Here the index v counts the spin (s) and the density 
{p) excitations. The phenomenological parameters Ky 
and completely characterize the low energy physics. 
For exactly solvable models like the ID lattice model, 
can be directly expressed in terms of the microscopic pa- 
rameters of the theory The gradients of the phases 
d(j)v are the density and the spin density fluctuations re- 
spectively. The canonical momenta H^ conjugated to the 
phases (j)^ are related to the spin and the density currents 

— ViiK^Wi,. In a rotationally invariant Fermi-gas the 
quantity Ks = 1, so that only independent parameters 
are Kp, and v^^p^ 0. 

As apparent from Eq. ([I]) a distinctive feature of the 
LM Hamiltonian is the complete separation of the spin 
and the "charge" degrees of freedom. In a spatially homo- 
geneous gas the spin and the charge waves propagate at 
the velocities and Vs respectively. This is the essence 
of the spin-charge separation phenomena. In our cold 
atoms system the charge and the spin waves are modeled 
by excitations of the total and the relative densities of 
the components. 

The local density approach to model a trapped gas as- 
sumes that the size of the atom cloud kp^, i.e. the 
size of the gas sample is much larger than the interparti- 
cle separation, consistent with ^ 1. The variation of 
and Vp, is assumed to originate only from the spatial 
dependence of the gas density, Kp[n] —^ Kp[n{x)] and 
Vp[n\ — > Vp[n{x)]. We will also assume that the numbers 



of particles of the both "spins" are the same (i.e. the 
total spin of the system is zero) . Our goal is now to de- 
rive the excitation spectrum in the charge and the spin 
sector analytically in the weak and strong coupling lim- 
its, and study the intermediate regime using numerical 
techniques. As a first step, this requires the calculation 
of the ground state density distribution. 

Within the local density approximation, the ground 
state of the system can be characterized using the 
Thomas-Fermi equilibrium condition: 
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Here E{n) is the internal energy of the gas per unit length 
of the e function of its (total) density, ^ is the 

chemical potential, V{x) = mtiPx'^ /2 is the longitudinal 
external potential, and w is the frequency of the longitu- 
dinal confinement. This equation is just the expression 
of the fact that the energy cost of adding a particle to the 
system equals to the chemical potential corrected by the 
local value of the external potential. Generally speaking, 
the external potentials acting on the two different " spin" 
components can be different. This feature can be used to 
generate offsets in the densities of the components and 
hence produce the spin and the density excitations with 
laser light (see below). 

In preparation for the interacting case we consider first 
a free gas ground state (g = 0). The density of the gas at 
a given position x is related to the local value of Fermi 
momentum, kp{x) — Tm{x)/2. The internal energy of 
the gas is just the density of the kinetic energy (the so 
called quantum pressure) E{n) — h'^7r^n^{x)/ 24771. Sub- 
stituting these expressions into Eq.(||) we find 
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for < Rtf, and otherwise. Here tiq = 7i{x = 0) = 
{Sfim/fT^Tr^y^^ is the density in the center and Rtf — 
(2/i/mw^)^/^ is the Thomas-Fermi size of the cloud . 
From the requirement that the integrated density equa s 
the particle number we have the condition /i — hLuN/2. 

The excitations of the gas can be found from the equa- 
tions of motion following from the Hamiltonian ([^): 



= K^{x)v^{x)Il„, Il„ 



d v^{x) d 
dx K^{x) dx 



(4) 



In particular, for an ideal gas we can use the density 
profile (^) and find the equations for the mode functions: 

where x = x/Rtf- The solution is given by 4>^n — 
A,^ sin(e„ arccosS) -|- i?^ cos(e„ arccosi). The discrete 
spectrum of eigenfrequencies is found by analyzing the 
boundary conditions: eu{n) = hio{7i-\- 1) both for the 
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spin and the density modes. The first modes (n = 0) 
and the wavefunctions (j)^ ~ cc/Vl — correspond to 
harmonic oscillations of the center of mass of the total 
density and the total spin (dipole modes). 

Before starting with the perturbation theory in small 
interaction parameter ^ ^ 1, we note that in a finite sys- 
tem there is additional energy scale, which is the level 
spacing [huj). In order for the interaction effects man- 
ifest themselves in a way similar to a bulk system, we 
need the interaction to be stronger than the level spac- 
ing, n{x)g 3> huj, in contrast to the case of very weak 
interactions studied in ||] . 

In a homogeneous gas the Luttinger parameters in the 
Hamiltonian (|^) to the lowest order in ^ = g/irhvp ^ 1 
can be found using perturbation theory: Kg = 1, Vg = 
vf{1 - ?/2), Kp = l- e/2 and Vp = vf{1 + ^/2) [||. 
The energy of the ground state can be obtained by aver- 
aging the interparticle interaction over the ground state, 
Eq = h'^n'^n'^ /24m + gn'^/A. In the spirit of the local 
density approximation we substitute a spatially depen- 
dent density n{x) in the expressions for homogeneous 
gas. Then, using Eq.(|^) we find that in the first or- 
der in ^ = g/TTVpiO), the density of the gas uniformly 
decreases by Sn{x) — —2gm/h^i:^ , i.e. the interaction 
reduces the density, as expected. This simple conclusion 
holds everywhere as long as g ^ h^n{x)/'m, i.e. 

[Rtf - x)/Rtf < [gm/h^nof ~ 0(^2) (5) 
The velocities of the spin and the density waves are 
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where As = ?> and Ap — 1. Using the expansion in pow- 
ers of ^ of the Luttinger parameters and the density pro- 
file (||) we find, that the frequency of the density dipole 
mode does not depend on the interaction (as it should 
be), while the the spin dipole mode shift is given by an 
integral logarithmically diverging at the border of the gas 
cloud. The divergence occurs due to localization of the 
excitations of a free gas close to the gas cloud border and 
arises in any potential, which is a power law in x. Using 
the condition (|[) to cut off the divergence, we find 



Sgm 



TT^h'^nTFiO) 



log 



Trh^riTFi^) 
2gm 



This shift is negative and can be observed by compar- 
ing the spin and the density oscillations of the gas cloud. 
Note that in a harmonic trap the perturbation theory re- 
quirement is stronger than in a homogeneous Luttinger 
liquid: we have to require ^ log(l/^) ^ 1 instead of sim- 
ply ? < 1- 

For higher modes the application of the perturbation 
theory in Eqs.(^ turns inconvenient and the frequencies 
of the excitations can be analyzed within the WKB ap- 
proximation. The accuracy of the WKB spectrum esti- 
mation is Xj-K^T? whereas the expected correc- 
tions are of order g/vF- Therefore for sufficiently high 
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Fig. 3 The level spacing (in units of uo) between the spin 
(solid line) and the "charge" (dashed line) modes vs. log ^ 
for different central densities: n(0) = 0.25 (solid), 
n(0) = 0.42 (dot-dashed), n(0) = 0.58 (dotted). 

n the eigenfrequencies can be reliably obtained from the 
WKB quantization condition p7|]. 



p{x)dx — h-K{n + a). 



(7) 



where p{x) is the WKB momentum corresponding to a 
given energy, n is the (integer) quantum number, xo is 
the classical turning point and the constant a = 1 is fixed 
by comparing the WKB results and the exact solutions 
of Eqs.(Q) for a weakly interacting gas. Substituting the 
dispersion relation e = v p^s{x)p{x) with the velocities (^) 
into Eq. (|^) we obtain the the same sort of logarithmically 
diverging integrals as those for in the perturbation theory 
above. By regularizing them using the condition (|^) we 
find, that 
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This simple WKB calculation confirms the interaction 
dependent split of the spin and the density oscillation 
frequencies. 

In the limit of very large interaction strength [g ^ 
■nhvF) the repulsion between the atoms of the two differ- 
ent species is very strong. Hence, the properties of the gas 
are similar to those of an ideal single component gas of in- 
distinguishable particles. The density profile is still given 
by Eq.(^), but now with no = Uoo = {2^im/h^'!r'^)^^'^ , 
fi = hNuj, and Rtf = i?oo = (2^/mu;2)i/2_ xhis dis- 
tribution is less dense and thus broader than that for a 
weakly interacting gas. The density wave speed is equal 
to the Fermi velocity vf — h-Krioo/rn and, after integra- 
tion in Eq.(|^), we find that spectrum of the density waves 
is the same as in the non-interacting case above. In turn, 
the relation between the energy and the WKB momen- 
tum for the spin wave is given by 
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where the coefRcient i? ~ 13 (±2) was found numeri- 
cally. Once again, using the quantization condition (|^) 
and cutting off the logarithmically divergent integral at 
the point n <^ g, we find, that 

ens — fuj->{n + ot) 



gm log{gm/ h^Uoc) ' 

where a ~ 1. As it is clear from the latter expression, 
the interaction profoundly changes the properties of the 
spin mode. In the limit of the strong interaction the level 
spacing decreases and is much smaller than that between 
the density waves (w). 

In order to confirm our analytical results, we performed 
a numerical calculation valid for arbitrary interaction 
strength based on a lattice model. Using the exact solu- 
tion IjlS , 19 for calculation of the Luttinger constants and 
the Thomas- Fermi approximation , we determined the 
WKB level spacings for the spin and the charge modes. 
The results are presented in Fig. 2 as a plot of the ex- 
citations level spacing vs. the dimensionless interaction 
strength g/nvp calculated at the center of the trap. 

As outlined in Fig. 1, wavepackets of the spin and den- 
sity excitations can be generated by short off-resonant 
and state selective laser pulses focused to a spot size £ 
with R ^ £ ^ kp^, where Rtf is the size of the atom 
cloud and kp^ the interparticle distance. This procedure 
is analogous to the MIT setup originally used to study 
propagation of sound waves in elongated condensates 
pO{ . Fig. lb shows as an example the wave packet dy- 
namics for the states \F — 1/2, Mp — ±1/2) of 6Li with 
interaction parameter ^ = 1, corresponding to N = 500 
particles at trap frequency of a; = IHz, = 250kIIz 
and scattering length as — 23A. Tuning near the Fesh- 
bach resonance (at B = 800G ) allows an increase of the 
scattering length by one order magnitude, allowing for 
N — 1000 atoms at a trap frequency of uj± = lOOkHz. 

Let us finally now turn to a discussion of the life time 
of the excitations and temperature requirements. The 
Hamiltonian (|^) represents only the first term in hydro- 
dynamic expansion q/kp <^ 1. The higher order terms 



originate from, for example, non-Hnearity of the fermionic 
spectrum and mix the excitations with each other. The 
first corrections are of the third power in 11 and d(f>, and 
hence lead to scattering of the excitations. To study the 
relaxation phenomena we switch for simplicity to a case 
of a single-component LM. According to Haldane, the 
generic term is given by Vint = jh'^ /m{d(j>)^ , where 7 '--^ 1 
Ig. The damping of the oscillations can be found from 
Fermi's Golden rule applied to the interaction Hamilto- 
nian Vint and using the second quantized representation 
for the phonon field. A straightforward calculation gives 
Ft^o ~ h?q*L/m^u. for fc^T ^ hutg, which is the decay 
rate in a process where a particle with the energy Hujq de- 
cays into a pair of particles with wi, 0^2 < L^q)- Note that 
cijg 3> w, otherwise there are no final states for such decay 
instability (i.e. the lowest excitations are very stable at 
very low temperatures). In the case of fc^T 3> hujqCbnt 
still T < Tp) we find Ft ~ q^Lk^T'^ /m^u^ which cor- 
responds to Landau damping, i.e. the contribution of 
a process in which the damping occurs by scattering a 
high-frequency excitation with hjj'^ ~ ksT. Both results 
contain the size of the sample L and are only valid in 
the collisionless regime F/w^ <C 1. This is the case for 
sufficiently small temperatures (or high number of par- 
ticles). Indeed, for the lowest excitations Lmu/h ^ N, 
ujq ^ uj ep/hN, so that T/uq ^ {ksT/epY < 1, 
i.e. the excitations are only weakly damped. 

In conclusion we performed the analysis of a double 
component Fermi gas confined in a harmonic trap. Based 
on the LM we have investigated the nature of the exci- 
tations and analyzed an experiment where spin-charge 
separation can be observed "directly" in experiments ad- 
dressing the spectral properties of the lowest excitations 
with laser light. 

Discussions with J.I. Cirac, J. von Delft, D. Jaksch 
and U. SchoUwock are gratefully acknowledged. Work 
supported in part by the Austrian Science Foundation 
and EU Networks, and the Institute for Quantum Infor- 
mation. 
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